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A SATELLITE-BASED LOCATION SYSTEM EMPLOYING 
DYNAMIC INTEGRATION TECHNIQUES 

RELATED APPLICATIONS 

Related subject natter is disclosed in the foUoAving applications and assigned to 
the same Assignee hereof: U.S. Patbit Application Serial No. 08/927,434 entitled "An Auxiliary 
System For Assisting A Wireless Temiinal In Determining Its Position From Signals Transmitted 
From A Navigation Satellite," inventors Robert Ellis Richton and Giovanni Vannucci; U.S. 
Patent Application Serial No.08/927,432\entitled "Telecommunications Assisted Satellite 
Positioning System," inventors Giovanni Vaimucci; U.S. Patent Application Serial No. 
09/321,075 entitled "Wireless Assisted GPS Using A Reference Location," inventors Robert 
Ellis Richton and Giovanni Vannucci; and U.SV Patent Application Serial No. 60/1 14,491 
entitled "Wireless Assisted Satellite Location Usiag a Reference Point," inventors Robert Ellis 
Richton and Giovanni Vannucci. Related subject matter is disclosed in the following application 
filed concurrently herewith and assigned to the same^ssignee hereof: U.S. patent application 
entitled "Satellite-Based Location System Employing Knowledge-Based Sequential Signal 
Search Strategy", Serial No. . \ 

FIELD OF THE INVENTION 

The present invention relates generally to wireless communication systems and, 
in particxxlar, to satellite-based location systems. 

BACKGROUND OF THE RELATED ART 

Satellite-based navigational systems provide accurate, three dimensional position 
information to worldwide users. Prior art satellite-based navigational systems, however, utilize a 
time consimiing search process for determining position information. Time consuming search 
processes are imdesirable in navigational systems particiilarly when the \iser is moving or in an 
emergency situation requiring immediate assistance. 

FIG. 1 depicts a well-known satellite-based navigational system referred to as 
Global Positioning System (GPS) 10. GPS 10 comprises a plurality of satellites l2-j and at least 
one GPS receiver 14, where/=l,2,. . .,n. Each satellite 12-/ orbiting earth at a known speed Vj and 
being a known distance apart from the other satellites l2-j. Each satellite 12-y transmits a GPS 
signal 1 l-y which includes a carrier signal with a known frequency / modulated xising a unique 
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pseudo-random noise (PN-y) code and navigational data (ND-y) associated with the particular 
satellite I2-7, wherein the PN-y code includes a unique sequence of FN chips and navigation data 
ND-y includes a satellite identifier, timing information and orbital data, such as elevation angle aj 
and azimuth angle ^j, FIG. 2 depicts a typical 20 ms frame of the GFS signal 1 1-y which 
5 comprises twenty fixll sequences of a FN-y code in addition to a sequence of navigation data 
ND-y. 

GPS receiver 14 comprises an antenna 15 for receiving GPS signals 1 1-y, a 
plurality of correlators 16-^ for detecting GFS signals 1 1-y and a processor 17 having software 
for determining a position using the navigation data ND-y, where ^ 1 ,2, . . . ,m. GPS receiver 1 4 

10 detects GPS signals 1 1-y via PN-y codes. Detecting GFS signals 12-y involves a correlation 
process wherein correlators l6-k are used to search for PN-y codes in a carrier frequency 
dimension and a code phase dimension. Such correlation process is inq)lemented as a real-time 
multiplication of a phase shifted replicated PN-y codes modulated onto a replicated carrier signal 
with the received GPS signals 1 1-y, followed by an integration and dmnp process. 

15 In the carrier frequency dimension, GPS receiver 14 replicates carrier signals to 

match the frequencies of the GFS signals 1 1-y as they arrive at GFS receiver 14. However, due to 
the Doppler effect, the frequency / at which GPS signals 1 1-y are transmitted changes an 
imknown amoimt before GPS signal 1 1-y arrives at GPS receiver 14 - that is, each GFS signal 
1 1-y should have a frequency when it arrives at GPS receiver 14. To account for the 

20 Doppler effect, GPS receiver 14 replicates the carrier signals across a frequency spectrum^^ 
ranging from /*-4/mm to /*-A^ until the frequency of the replicated carrier signal matches the 
frequency of the received GPS signal 1 1-y, wherein 4/min and are a minimxmi and maximum 
change in frequency GFS signals 1 1-y will undergo due to the Doppler effect as they travel from 
satellites 12-y to GPS receiver 14, i.e., 4/mm ^ 45 ^ A/max- 

25 In the code phase dimension, GFS receiver 14 replicates the unique PN-y codes 

associated with each satellite 12-y. The phases of the repUcated PN-y codes are shifted across 
code phase spectrums R^<spec) until replicated carrier signals modulated with the replicated PN-y 
codes correlate, if at all, with GFS signals 1 1-y being received by GPS receiver 14, wherein each 
code phase spectrum R,{spec) includes every possible phase shift for the associated FN-y code. 

30 When GFS signals 1 1-y are detected by correlators 16-^, GPS receiver 14 extracts the navigation 
data ND-y from the detected GPS signals 1 1-y and uses the navigation data ND-y to determine a 
location for GPS receiver 14, as is well-known in the art. 
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Correlators 16-^ are configured to perform parallel searches for a plurality of 
PN-y codes across the frequency spectrum^^pcc and the code phase spectnims R,{spec). In other 
words, each of the plurality of correlators 16-^ are dedicated to searching for a paiticxilar PN-j 
code across each possible frequency between />-4/mm to /t-A/max and each possible for that PN-y 
5 code. When a correlator IS-k conq)letes its search for a PN-y code, the correlator 16-^ searches 
for another PN-y code across each possible frequency between /^-A/mm to /J-A/max and each 
possible phase shift for that PN-y code. This process continues until all PN-y codes are 
collectively searched for by the plurality of correlators IS-k, For exan^le, suppose there are 
twelve satellites 12-y, thus there would be twelve unique PN-y codes. If GPS receiver 14 has six 

10 correlators 16-/:, then GPS receiver 14 would use its correlators 16-^ to search for two sets of six 
different PN-y codes at a time. Specifically, correlators l6-k search for the first six PN-y codes, 
i.e., correlator 16-1 searches for PN-1, correlator 16-2 searches for PN-2, etc. Upon convicting 
the search for the first six PN-y codes, correlators l6-k search for the next six PN-y codes, i.e., 
correlator 16-1 searches for PN-7, correlator 16-2 searches for PN-8, etc. 

15 For each PN-y code being searched, correlator 16-A: performs an integration and 

dump process for each combination of frequency and phase shifts for that PN-y code. For 
example, suppose the frequency spectrum,^ includes 50 possible frequencies for the carrier 
signal and the code phase spectrum R,<spec) for a PN-y code includes 2,046 possible half-chip 
phase shifts. To search for every possible combination of frequency and half-chip phase shifts 

20 for the PN-y code, the correlator 16-^ would then need to perform 102,300 integrations. A typical 
integration time for correlators l6-k is 1 ms, which is generally sufficient for GPS receiver 14 to 
detect GPS signals 1 1-y when antenna 15 has a clear view of the sky or a direct line-of-sight to 
satellites 12-y. Thus, for the above example, 102 .3 seconds would be required for one correlator 
16-Jt to search every possible combination of frequency and half-chip phase shifts for a PN-y 

25 code. 

GPS receivers, however, are now being incorporated into mobile-telephones or 
other types of mobile communication devices which do not always have a clear view of the sky. 
Thus, GPS receiver 14 will not always have a clear view of the sky. In this situation, the 
signal-to-noise ratios of GPS signals 1 1-y received by GPS receiver 14 are typically much lower 
30 than when GPS receiver 14 does have a clear view of the sky, thus making it more difficult for 
GPS receiver 14 to detect the GPS signals 1 1-y. To compensate for weaker signal-to-noise ratios 
and enhance detection of GPS signals 1 1-y, correlators I6'k can be configured with longer 
integration times. In other words, detection is more accxirate with longer integration times. A 
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sufficient integration time, in this case, would be approximately 1 second. Thus, for the exanq)le 
above, 102,300 seconds would be required for a correlator 16-A: to search for every possible 
combination of frequency and half-chip phase shifts for a PN-y code. Longer integration times 
result in longer acquisition times for detecting GPS signals 1 l-y. Longer acquisition times are 
5 undesirable. 

Wireless assisted GPS (WAG) systems were developed to facilitate detection of 
GPS signals 1 l-j by GPS receivers configured with short or long integration times. The WAG 
system facilitates detection of GPS signals 1 l-y by reducing the number of integrations to be 
performed by correlators searching for GPS signals 1 1-j. The nimiber of integrations is reduced 

10 by narrowing the frequency range and code phase ranges to be searched. Specifically, the WAG 
system limits the search for GPS signals 1 l-y to a specific frequency or frequencies and to a 
range of code phases less than the code phase spectrum R,<spec). 

FIG. 3 depicts a WAG system 20 con:q)rising a WAG server 22, a plurality of 
base stations 23 and at least one WAG client 24. WAG server 22 includes a GPS receiver 26 

15 having an antenna 27 installed in a known stationary location with a clear view of the sky. GPS 
receiver 26 would typically have correlators configured with short integration times because 
antenna 27 has a clear view of the sky. WAG server 22 being operable to communicate with 
base stations 23 either via a wired or wireless interface. Each base station 23 has a known 
location and provides communication services to WAG cUents located within a geographical area 

20 or cell 25 associated with the base station 23, wherein each cell 25 is a known size and is divided 
into a plurality of sectors. WAG client 24 includes a GPS receiver 28 and perhaps a 
mobile-telephone 27, and is typically in motion and/or in an unknown location with or without a 
clear view of the sky. GPS receiver 28 having correlators typically configured with long 
integration times. Note that the term "mobile-telephone," for purposes of this application, shall 

25 be construed to include, but is not limited to, any communication device. 

FIG. 4 is a flowchart 300 illustrating the operation of WAG system 20. In step 
310, WAG server 22 detects a plurahty of satellites 12-y via their GPS signals 1 \-j using its GPS 
receiver 26. WAG server 22 acquires the following information from each detected satellite 
12-7": the identity of satellite 12-y and frequency J}, code phase, elevation angle q, and azimuth 

30 angle (J), associated with the detected satellite 12-7, wherein the elevation angle a,- is defined as 
the angle between the line of sight from WAG server 22 or client 24 to a satellite 12-y and a 
projection of the line of sight on the horizontal plane, and the azimuth angle is defined as the 
angle between the projection of the line of sight on the horizontal plane and a projection of the 
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north direction on the horizontal plane. See FIG. 5, which depicts an elevation angle q, and an 
azimuth angle <t)y corresponding to a satellite l2-j and a WAG server 22 or WAG client 24. 

In step 315, WAG server 22 receives sector information from base station 23 
currently in communication with or serving WAG client 24, wherein the sector information 
5 indicates a sector WAG client 24 is currently located. In step 320, WAG server 22 makes an 
initial estimate of WAG client's position based on the known location of the serving base station, 
the cell size associated with the serving base station, and the sector in which WAG client 24 is 
currently located. In one embodiment, WAG server 22 initially estimates that WAG client 24 is 
located at a reference point within the sector, e.g., point at approximate center of sector. In 

10 another embodiment, WAG server 22 initially estimates WAG client 24 's position using 
well-known forward link triangulation techniques. 

In step 330, for each detected satellite 12-y, WAG server 22 uses the information 
acquired from the detected GPS signals 1 l-j to predict a frequency j^{r) at the reference point and 
a code phase search range R^{sect) which includes all possible code phases for GPS signal ll-j 

15 arriving anywhere within the sector where WAG client 24 is currently located. In step 340, 
WAG server 22 transmits a search message to the serving base station 23, wherein the search 
message includes, for each detected satellite I2'j\ information regarding the associated PN-y 
code, predicted frequency Jj{t) and code phase search range Resect). 

In step 350, serving base station 23 transmits the search message to WAG client 

20 24 which, in step 360, begins a parallel search for the satellites ll-j indicated in the search 

message. Specifically, WAG client 24 will use its correlators to simultaneoxisly search for each 
of the GPS signals 1 1-y at the predicted frequency j^{r) within the limitations of the code phase 
search range R^<sect) indicated in the search message. Thus, the number of integrations is 
reduced to the predicted frequency j^<r) within the limitations of the code phase search range 

25 R/(sect). However, due to the long integration times of GPS receiver 28 in WAG client 24, the 
search time is still considered time consuming. Accordingly, there exists a need to facilitate 
detection of satellites I2-7. 

SUMMARY OF THE INVENTION 
30 The present invention is a method and apparatus for facilitating detection of 

satellite signals using a dynamic integration technique. The present invention uses a dynamic 
integration technique in which integration time periods of correlators are adjusted according to 
signal strength measxirements of satellite signals received at GPS receivers. Specifically, 
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integration time periods are inversely adjxisted, either proportionally or non-proportionally, to 
received strengths of signals being searched. In one embodiment, when received signal strengths 
are above or equal to a threshold value, the integration time period is short, e.g., 1 millisecond. 
By contrast, when received signal strengths are below the threshold value, the integration time 
5 period is long, e.g., 1 second. For exanaple, when WAG cUent is indoors and the satellite signals 
are weak, the integration time period is 1 second, and when WAG cUent is outdoors and the 
satelUte signals are strong, the integration time period is 1 millisecond. Advantageously, because 
the integration time periods are dynamically and inversely adjusted according to received signal 
strengths of satellite signals, detection of satellite signals is facilitated without sacrificing 
10 accuracy of detection. To fiirther facilitate detection, one embodiment uses a sequential search 
technique during which satelUtes are searched for in a sequential marmer based on knowledge 
acquired fi-om detecting other satellite signals. 

BRIEF DESCRIPTION OF THE DRAWINGS 
15 The features, aspects, and advantages of the present invention will become better 

understood with regard to the following description, appended claims, and accompanying 
drawings where 

FIG. 1 depicts a well-known satellite-based navigational system referred to as Global 
Positioning System (GPS); 
20 FIG. 2 depicts a typical 20 ms firame of a GPS signal; 

FIG. 3 depicts a Wireless Assisted GPS (WAG) system; 

FIG. 4 depicts a flowchart illustrating the operation of the WAG system of FIG. 3; 
FIG. 5 depicts an elevation angle a, and an azimuth angle corresponding to a satellite 
and a WAG server or WAG client; 
25 FIG. 6 depicts a flowchart illustrating a sequential search technique and dynamic 

integration technique used in accordance with one embodiment of the present invention; 

FIG. 7 depicts a chart illustrating two curves that may be used for determining 
integration time periods; and 

FIG. 8 depicts a sanqjle power spectrum density of auto correlation for eight correlators 
30 spaced half a chip period apart. 
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DETAILED DESCRIPTION 

The present invention is described herein with reference to the aforementioned 
WAG system. This should not be construed, however, to limit the present invention in any 
manner. 

5 FIG. 6 is a flowchart 600 illustrating a sequential search technique used in 

accordance with one embodiment of the present invention. In step 605, WAG client 24 receives 
a search message from its serving base station 23 or WAG server 22. In the present invention, 
the search message includes, for each satellite ll-j detected by WAG server 22, information 
regarding the associated FN-j code, predicted frequency ^{r) at a reference point within the 
10 sector/cell where WAG client 24 is currently located, code phase search range Resect) including 
all possible phase shifts for a GPS signal 1 l-J transmitted by satellite 12-y and arriving within the 
sector/cell where WAG client 24 is currently located, and orbital data including elevation angle 
aj and azimuth angle ^j. 

In step 610, WAG client 24 selects a first satellite I2~j indicated in the search 
W 15 message to search. WAG client 24 uses one or more criteria in a set of first satellite selection 
U criteria to select the first satellite 12-y. In one embodiment, the set of first satellite selection 

;r! criteria are as follows: (1) maximize utilization of correlators; (2) minimize search time; and (3) 

IjJ maximize the amount of information regarding location of WAG client 24 (or GPS receiver 14 or 

antenna 15). The first criteria of maximizing utilization of correlators involves using as many of 
=0 20 the available correlators to simultaneously search for a satellite l2-j\ The second criteria of 
;rl minimizing search time involves reducing the number of integrations to be performed by each 

^fl correlator, e.g., each correlator performs one integration. Reducing the number of integrations to 

be performed by each correlator essentially means selecting a satellite 12-y having the smallest 
associated code phase search range R/ indicated in the search message. 
25 The third criteria of maximizing the amoimt of information regarding the 

location of WAG client 24 involves selecting a satellite 12-y that, when detected, will indicate an 
area in the sector where WAG client 24 is located. For exanq)le, a satellite l2-j with a small 
elevation angle a^, when detected, will indicate a narrow strait in the sector where WAG client is 
located, whereas a satellite l2-j with a large elevation angle a,- will indicate a wider strait in the 
30 sector where WAG client is located. 

Upon selecting a first sateUite 12^' to be searched, in step 620, WAG cUent 24 
determines an integration time period for correlators 16-^ at which to search for the first sateUite 
l2-j, wherein the integration time period will depend upon the received signal strength 
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measurements at WAG client 24 (or server 22) of signals lUj transmitted by the first or any 
satellite 12-y. 

In one embodiment, WAG client 24 measxires the signal strength of signals 
transmitted on or about fi-equency f^^j as indicated in the search message for the first sateUite 
12-y. In another embodiment, WAG client measures the signal strength of signals transmitted on 
firequency / 

The integration time periods vary inversely, either proportionally or non- 
proportionally, to the signal strength measurements of signals 1 l-j. For exanq>le, if the received 
signal strength measurements of signals I l-j increase, the integration time periods for correlators 
l6-k decrease. In one embodiment, if the received signal strength measxirements of signals 1 1-^ 
at WAG client 24 is below a threshold value, the integration time period is maximized, and vice- 
versa. For example, if WAG client 24 measures the strength of signals 1 l-j to be below the 
threshold value, the integration time period of correlators 16-^ is set to I second, whereas if 
WAG client 24 measures the strength of signals 1 l-j to be above or equal to the threshold value, 
the integration time period of correlators 16-^ is set to 1 nfiillisecond. 

In another embodiment, there are multiple threshold values and integration time 
periods for correlators 16-^. Suppose there are two threshold values and three durations of 
integration time periods. If the signal strength is below the lowest threshold value, the 
integration time period is set to shortest duration; if the signal strength is above the lowest 
threshold value but below the highest threshold value, the integration time period is set to second 
shortest or second longest duration; and if the signal strength is above the highest threshold 
value, the integration time period is set to longest duration. In other embodiments, the 
integration time periods are determined according to the received signal strengths and a curve or 
mathematical equation. See FIG. 7, which depicts a chart 70 illustrating two possible 
relationships between received signal strength measurements and integration time periods. 

In yet another embodiment, relative signal strengths are used to determine 
whether GPS receiver 14 is indoors or outdoors which, in turn, determines the integration time 
periods. For example, suppose GPS receiver 14 has eight correlators which are spaced half a 
chip period apart. Each correlators samples or measures signal strengths (for a particular signal 
1 l-j or PN-y sequence) every millisecond over a twenty millisecond time period, i.e., each 
correlator perform twenty samples. For each set of san:5)les collected by a correlator, a Fourier 
transform is applied to obtain a power spectrum density of auto correlation for each correlator. 
FIG. 8 depicts a three dimensional chart 80 illustrating a possible power spectrum density of auto 
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correlation for eight correlators spaced half a chip period apart. A maximum power spectrum 
density value over all the correlators (Pmax) is divided over an average power spectrum density 
value (Pavg) over all the correlators to obtain a power spectrum density ratio (Pratio). 

J^max 



fO.-^ Pavg 



= Pr, 



5 The power spectrum density ratio Piatio is then conq)ared to a threshold value to 

determined integration time periods. Basically, a small Pratio would indicate that GPS receiver 14 
is indoors, whereas a large or larger Pratio would indicate that GPS receiver 14 is outdoors. 
Specifically, when GPS receiver 14 is indoors, more multi-path signals are present resulting in a 
wider distribution of a signal's energy. Thus, Pmax should not be much greater than Pavgand Pratio 

10 should be a small value (e.g., a value greater than but close to one). By contrast, when GPS 

receiver 14 is outdoors, less multi-path signals are present resulting in more of a concentration of 
a signal's energy onto a single signal. Thus, Pj^x should be much greater than Pavg and Pratio should 
be a large or larger value (e.g., a value greater than and not close to one). 

In step 625, WAG client 24 searches for the first satellite ll-j for the determined 

1 5 integration time period using the frequency j^<r) and code phase search range R^<sect) indicated in 
the search message for the first satellite 12-7*. Once the first satellite 12-y has been detected, in 
step 630, WAG client 24 predicts a first area in which WAG cUent 24 may be located using 
information extracted from a GPS signal 1 l-j transmitted by the first satellite l2-j\ as is well- 
known in the art. The first predicted area typically being a strait or small area within the sector 

20 where WAG client 24 is currently located. Such calculation is later used to narrow down the 
code phase search range of subsequent satellite searches. 

In step 640, WAG cUent 24 uses the search message to pick a second satellite 
12-y to search. WAG client 24 uses one or more criteria in a set of second sateUite selection 
criteria to select the second satellite l2-j. In one embodiment, the set of second satellite selection 

25 criteria are as follows: (1) maximize utilization of correlators; (2) minimize search time; and (3) 
maximize the amoimt of additional information regarding location of WAG client 24 (or GPS 
receiver 14 or antenna 15) when used in conjunction with the resxilts of the first search. The first 
and second criteria being identical to the first and second criteria of step 610. The third criteria 
involves selecting a second satellite 12-y that will result in an area which intersects least, but 

30 nevertheless intersects, with the first predicted area. In one embodiment, the second satellite 12-y 
selected is a satellite 12-y that forms an angle of approximately 90° with the first satellite and 
WAG server 22 or WAG client 24, wherein WAG server 22 or WAG chent 24 is the vertex. The 
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angle between the fist and second satellites and WAG server 22 or client 24 can be determined 
using a difference between azimuth angles associated with the first and second satellites. 

In step 645, WAG client 24 redefines or narrows down the code phase search 
range R^{sect) indicated in the search message for the second satellite ll-j based on the first 
5 predicted area. The redefined or narrowed down code phase search range R^{sect) is hereinafter 
referred to as a '^predicted code phase search range R^pred)." The predicted code phase search 
range R,{pred) for the second satellite includes all possible phase shifts for GPS signals 1 l-j 
transmitted by the second satellite and arriving in the first predicted area. Since the first 
predicted area is a strait or small area within the sector where WAG client 24 is currently 
10 located, the corresponding predicted code phase search range R/(pred) will be narrower than the 
corresponding code phase search range R^{sect) originally indicated in the search message for the 
second satellite. 

In step 648, WAG client 24 determines an integration time period for correlators 
16-Jt at which to search for the second satellite 12-y based on received signal strength 

15 measurements. As in step 620, WAG client 24 can either measure the signal strength of signals 
transmitted on or about frequency f+^j as indicated in the search message for the second satellite 
12-y or on frequency / In step 650, WAG cUent 24 searches for the second satellite 12-y, for the 
integration time period determined in step 638, using the frequency j^{r) indicated in the search 
message and the predicted code phase search range R/(pred) for the second satellite 12-y. In 

20 another embodiment, the integration time period may be re-determined prior to step 650. Once 
the second satellite 12-y has been detected, in step 655, WAG client 24 predicts a second area in 
which WAG client 24 may be located xising information extracted from a GPS signal 1 1-y 
transmitted by the second satellite 12-y. Like the first predicted area, the second predicted area is 
typically a strait or small area within the sector in which WAG cUent 24 is currently located. 

25 The intersection of the first and second predicted areas effectively reduces the 

size of the search area in which WAG client 24 may be located. In step 660, WAG cUent 24 uses 
the intersected area and the code phase search range R^<sect) indicated in the search message to 
predict code phase search range R^<pred) for the remaining satellites 12-y indicated in the search 
message, thereby facilitating detection of the remaining satellites 12-y. Such predicted code 

30 phase search ranges R/(pred) include code phases for GPS signals 1 1-y transmitted by the 

remaining satellites indicated in the search message and arriving anywhere within the intersected 
area. 
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In step 665, WAG client 24 determines integration time periods for correlators 
16-^ at which to search for the remaining satellite 12-7 all other satellites 12-y indicated in 
the search message other than the first and second satellites) based on received signal strength 
measurements.. As in step 620, WAG client 24 can either measure the signal strength of signals 

5 transmitted on or about frequency f^d/j as indicated in the search message for the remaining 

satellites ll-j or on frequency / In step 670, WAG cUent 24 searches for the remaining satellites 
ll'j, for the integration time period determined in step 665, within the confines of the predicted 
code phase search ranges R/(pred) for the remaining satellites ll-j. In an embodiment of the 
present invention, WAG client 24 uses its correlators to perform parallel searches for two or 

10 more remaining satellites \2-j. Upon detecting the remaining satellites 12-7, step 680, WAG 
cUent 24 calculates its location using the navigation data ND-y extracted from GPS signals 1 1-y 
transmitted by at least three satellites 1 l-j, as is well-known in the art. 

The present invention is described herein with reference to certain embodiments, 
including an embodiment in which the first, second and all or some of the remaining satellites are 

15 searched sequentially. Other embodiments are possible. For example, the sequential search of 
the present invention may involve GPS receiver 28 searching in parallel for the first and second 
satellites, and then searching in parallel for all or some of the remaining satellites. The present 
invention is also applicable to non-GPS satellite-based or non-satellite-based navigation system. 
Accordingly, the present invention should not be linwted to the embodiments disclosed herein. 
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